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ABSTRACT 

Analysis of several recent ROSAT HRI observations of the gravitationally 
lensed system Q0957+561 has led to the detection at the 3a level of the cluster 
lens containing the primary galaxy Gl. The total mass was estimated by 
applying the equation of hydrostatic equilibrium to the detected hot intracluster 
gas for a range of cluster core radii, cluster sizes and for different values of the 
Hubble constant. 

X-ray estimates of the lensing cluster mass provide a means to determine 
the cluster contribution to the deflection of rays originating from the quasar 
Q0957+561. The present mass estimates were used to evaluate the convergence 
parameter k, the ratio of the local surface mass density of the cluster to 
the critical surface mass density for lensing. The convergence parameter, n, 
calculated in the vicinity of the lensed images, was found to range between 0.07 
and 0.21, depending on the assumed cluster core radius and cluster extent. This 
range of uncertainty in k does not include possible systematic errors arising 
from the estimation of the cluster temperature through the use of the cluster 
luminosity-temperature relation and the assumption of spherical symmetry of 
the cluster gas. Applying this range of values of k to the lensing model of 
Grogin & Narayan (1996) for Q0957+561 but not accounting for uncertainties 
in that model yields a range of values for the Hubble constant: 67 < H < 82 
km s -1 Mpc -1 , for a time delay of 1.1 years. 

Subject headings: gravitational lensing — quasars: individual (Q0957+561) — X- 
rays: galaxies 
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1. INTRODUCTION 

Refsdal (1964a 1964b 1966) outlined a method for determining a global value for 
the Hubble constant by measuring the light travel delay between different images of a 
gravitationally lensed system. Together with a model that describes the gravitational 
potential of the lens, the time delay defines the geometrical dimensions of the lens system 
and therefore may lead to the determination of the Hubble constant. 

The criteria that make a gravitational lens system (GLS) suitable for measuring the 
Hubble constant can be summarized as follows: (a) The sky positions and redshifts of 
the lens and of the images of the lensed object should be measurable accurately; (b) The 
mass distribution of the lens should be subject to accurate and reliable estimation. When 
the lens has multiple components, as for Q0957+561, this estimation becomes especially 
difficult and benefits from many different types of observations and also from many and 
extended images, such as rings. For example, if the types of observations include radio and 
optical emission lines, which are less sensitive to microlensing, flux ratio measurements are 
more useful in constraining the lens model; and (c) The differences in propagation time 
from the lensed object to us through (at least two) different images ("time delay") should 
be measurable with high accuracy. The source, i.e., the lensed object, must therefore be 
variable on time scales far shorter than the time delays, to allow accurate estimates to be 
made of the latter. 

Several potential candidates for determining H are the quasar Q0957+561 with two 
images and a measured time delay of 1.1 years (Pelt et al. 1996; Haarsma et al. 1997; 
Kundic et al. 1997), the quadruple image system PG 1115+080 with a measured time delay 
between components B and C of about 25 days (Schechter et al. 1997; Bar-Kana 1997), 
and B0218+357 with a time delay of 12 ± 3 days (Geiger & Schneider 1996; Corbett et al. 
1996). Combining measured time delays with a detailed lensing model, Refsdal's suggestion 
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has been applied to Q0957+561 (Falco, Gorenstein, & Shapiro 1991; Grogin & Narayan 
1996, hearafter GN) and PG1115+080 (Keeton & Kochanek 1997). 

In this paper we focus on a study of the properties of Q0957+561 (Walsh, Carswell 
& Weyman 1979), one of the most extensively studied lensing systems which has been 
monitored since its discovery in both the optical and radio (see, for example, Schild & 
Thomson 1995; Vanderriest et al. 1989; Haarsma et al. 1997). The primary lensing 
components are an elliptical cD galaxy at a redshift of 0.355, usually referred to as the Gl 
galaxy, a cluster of galaxies containing the Gl galaxy, and a group of galaxies at a redshift 
of 0.5. The most thorough optical spectroscopic and photometric study to date of the 
field around Q0957+561 has been presented by Angonin-Willaime, Soucail, & Vanderriest 
(1994). 

One of the major concerns in modeling Q0957+561 is a mass sheet degeneracy which 
arises from the presence of the cluster of galaxies at z=0.36. In particular if one were to 
modify the assumed radial mass density profile k{6) of the lens to \k{6) + (1 - A), where 1 
- A represents a constant surface mass density sheet term, then the new mass distribution 
would also satisfy the observational constraints. However the resulting value for H would 
be scaled by the constant 1 - A. Since A must always be positive one may always determine 
an upper limit on the Hubble constant for A = (Falco et al. 1991). 

As also noted by Falco et al. (1991) and later GN, this degeneracy may be broken if 
one were to measure the velocity dispersion or total mass distribution of either the cluster 
of galaxies or the principal lensing galaxy Gl. Recently Falco et al. (1997) have obtained a 
velocity dispersion for the central lens galaxy Gl of a v = 266 ± 12 km s" 1 or a v = 279 ± 12 
km s" 1 depending on the interpretation of the spectroscopic slit data. 

In this paper we focus on improving the cluster model and we present an estimate of 
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the total mass of this lensing cluster of galaxies through X-ray measurements of thermal 
emission from the intracluster gas. In section 2 we describe the spatial analysis of the 
X-ray data of Q0957+561 obtained from several ROSAT HRI observations and we present 
an estimate for the convergence parameter k(9) of the cluster. Section 3 describes the 
spectroscopic analysis of ROSAT position sensitive proportional counter (PSPC) and 
ASCA Gas Imaging Spectrometer (GIS) X-ray observations at different epochs in order to 
investigate the variability of different spectral components, and presents various scenarios to 
explain the observed variability in the X-ray flux. Finally section 4 presents our conclusions 
from the analysis of the X-ray observations of the Q0957+561 GL system. In particular, we 
discuss how our results allow bounds to be placed on the Hubble constant, albeit with these 
bounds still subject to systematic errors that we cannot yet evaluate accurately. 



2. SPATIAL ANALYSIS: A 3a DETECTION OF THE CLUSTER THAT 

CONTAINS THE Gl GALAXY 

The limits on the mass and luminosity of the cluster of galaxies that contains the 
primary lensing cD elliptical galaxy Gl presented in Chartas et al. (1995) were based on 
the ROSAT HRI observation of Q0957+561 from 1992 May 1. We now estimate the mass 
and luminosity of this cluster utilizing additional data obtained from the three subsequent 
ROSAT HRI observations in 1992 October 20, 1992 October 22, and 1994 April 11 (see 
Table 1). The October 20 ROSAT HRI observation is centered on the galaxy NGC 3079 
with the quasar Q0957+561 located 12 arcmin off-axis. Because of the dependence of the 
properties of the point spread function (PSF) with off-axis angle, we performed the analysis 
for this off-axis observation separately and combined the results with those obtained from 
the on-axis HRI observations of Q0957+561. To look for extended cluster emission centered 
on the galaxy Gl, we first merged the on-axis observations of 1994 April, 1992 October, 
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and 1991 May. The image files containing the data for these observations were shifted 
appropriately such that the center of Q0957+561 image B from each observation was 
co-aligned. The shifts were relatively small, of the order of a few arcsecs. Figure 1 shows 
the resulting merged image of the three HRI observations. X-ray emission surrounding the 
lensed images can be seen. To determine the significance of this extended emission and 
whether it originates from hot cluster gas emission or from the PSF wings of the lensed 
images, we searched for a statistically significant excess in detected counts originating from 
cluster emission around the Gl galaxy. In particular we determined the number of counts 
within annuli of inner radii ranging from 15" to 100" in steps of 5" and outer radius of 120" 
centered on Q0957+561 for the merged ROSAT HRI observations and separately for the 
off-axis observation. The background was taken from an annulus of inner and outer radii 
of 340" and 460". For inner radii less than about 25" the contribution from the PSF's of 
images A and B of Q0957+561 is significant whereas for inner radii greater than about 
60" the signal to noise drops considerably leading to a large uncertainty in the determined 
source counts, as seen in Figure 2. We modeled two point sources at quasar image locations 
A and B. Each modeled source PSF was normalized such that the number of counts within 
a radius of 3" from the center of the modeled source was equal to the number of counts 
within a radius of 3" of the observed source. The net counts N d (ri 2 ) detected by the 
ROSAT HRI originating from cluster emission within an annulus of inner radius r 1 and 
outer radius r 2 were determined by subtracting background and point source contributions 
and is given by the expression, 

N cl (r 12 ) = N tot (r 12 ) - N bkg (r 12 ) - N A (r 12 ) - N B (r 12 ) (1) 

where Nbkg{j\ 2 ) are the background counts in the annular region r± 2 and iV^r^), Nb{j\ 2 ) 
are the counts within the region r 12 originating from the A and B quasar images of 
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Q0957+561, respectively. 

To compute the total cluster counts, we estimated a correction term for the cluster 
counts beyond the range of the inner and outer radii of the annulus (as described in Chartas 
et al. 1995). The correction term depends on the adopted value for the cluster core radius 
and cluster radius. For our calculations we consider core radii ranging between r c = 5" and 
r c = 45" and cluster radii ranging between 120" and 280". For a typical cluster core radius 
of r c = 15" (0.067 Mpc, assuming H = 75 km s _1 Mpc -1 and q = 0), we find that the 
number of counts within a 20" circle centered on Q0957+561 is about 60% of the counts 
found in the annulus with inner radius of 20" and outer radius of 120", whereas for r c = 35" 
we find that about 18 % of the X-ray counts originate within 20". No X-ray sources 
detected within the HRI field of view lie within the annuli used as source and background 
extraction regions. For the four HRI observations of Q0957+561 with a total exposure time 
of 85154sec, we estimate that a total of about 300 ± 100 counts (with the quoted error 
being at the la level) originated from cluster emission and approximately 2900 counts from 
the quasar images. 

To convert the measured HRI counts into luminosities we have assumed a Raymond- 
Smith spectrum characterized by a temperature T e , a galactic column density Nh = 
0.87 x 10 20 cm -2 (Stark et al. 1992), 30% cosmic abundances (see, e.g., Henriksen 1985; 
Hughes et al. 1988; and Arnaud et al. 1987) and the L x - T e correlation determined by 
David et al. (1993). The relation between cluster luminosity and temperature for a value of 
the Hubble constant of H = lOOh km s" 1 Mpc" 1 is, T e (keV) = l(}- 12m54 (h 2 L x ) b with b = 
0.29. 

There is considerable scatter in the observational data used by David et al. (1993) 
to produce the L x vs. T e correlation. A recent study however by Markevitch (1998) has 
shown that the scatter in the L x - T e relation is greatly reduced when the emission from 
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cooling flow regions is excluded in the analysis. To derive the cluster mass, we assume a 
typical density profile p{r) of the hot gas (e.g., Jones & Forman 1996) given by : 



and apply the equation of hydrostatic equilibrium to the intracluster hot gas, which we 
assume to be isothermal. The total mass of the lens within a radius r is given by 



where r c is the core radius of the cluster, \im v is the mean molecular weight of the gas, and 
f3 is the ratio of kinetic energy per unit mass in galaxies to kinetic energy per unit mass in 
the gas. 

The derived cluster temperature depends on the assumed value of H and the assumed 
core radius of the cluster. Therefore, we determined the cluster mass within a radius of 
IMpc for a range of cluster core radii and for H values of 50 and 75 km s -1 Mpc" 1 . The 
results of this analysis are shown in Figure 3. 

We use the estimate of the cluster mass to evaluate the convergence parameter k, the 
ratio of the local surface mass density of the cluster to the critical surface mass density for 
lensing. The convergence parameter together with the time delay Atba of the GLS are two 
remaining principal parameters with the largest uncertainties that enter in the gravitational 
models recently developed by GN used to determine the Hubble constant. GN obtain, 




(2) 




(3) 



H = (85±f)(l - K)(l.lyr/AT B A)km S - 1 Mpc- 1 . 



(4) 



Recent deep HST observations of Q0957+561 by Bernstein et al. (1997) have detected 
several new features and provide improved positions for the lens components. In particular 
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the position of the center of Gl is found to lie at the position previously determined by 
VLBI (see, e.g., Falco et al. 1991) and significantly distant from the position used in the 
GN derivation of equation 4. The previously published estimates for k are very uncertain 
because they depend (strongly) on the assumption that the cluster potential may be 
represented as a softened isothermal sphere. The velocity dispersion and its relation to 
the mass of the cluster, the cluster core radius, the cluster shape and the position of the 
cluster center with respect to the lensed object also need to be known. The present HRI 
observations provide a partly independent, but still mo del- dependent, estimate of the total 
mass of the lens and hence of the convergence parameter k through the expression 



K(X) = ^ (5) 



where is the surface mass density of the cluster lens and S cr is the critical surface 
mass density (see, e.g., Blandford & Narayan 1992), the density the cluster must have to 
produce multiple images if it were the only lensing agent, and is given by: 

c 2 D os 

Scr = AnGD ds D od (6) 
where D os , D ds , and D od are the angular-diameter distances between source and observer, 
deflector and source, and deflector and observer, respectively. We evaluated the critical 
surface mass density for the gravitationally lensed system Q0957+561 to be S cr = 
0.9452(0.9224)/i gr cm" 2 , for q = 0(0.5). 

In estimating H Q from the GLS, we need to know the cluster contribution to the lensing 
of the quasar. In this treatment, we linearly superpose the effects of Gl and the cluster. 
We first view Gl's contribution to the lensing as producing a pair of images, A' and B', as 
shown in Figures 4 and 5. The contribution of the central galaxy Gl has been represented 
by SPLS and FGS models (GN, Falco et al. 1991). 
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For the cluster, we assume a spherically symmetric potential which produces an 
outward radial offset of the lensed components with respect to the cluster center (see 
Figures 4 and 5). To simplify the calculations we incorporate the thin lens approximation 
in which the cluster is projected onto a plane perpendicular to the line of sight, on the 
plane of the sky. 

To calculate the angle, a, by which a ray is deflected by the cluster lens, we integrated 
the deflection over the entire cluster. The scaled deflection angle (Schneider, Ehlers, & 
Falco 1992) is given by 



i -* -*r 

1 ' j2 / , ^ X — X 



a(x) = - d 2 x'K(x>) (7) 

7T J \X — X Y 

where x' represents the scaled location of a projected mass element in the lens plane with 
respect to the cluster center and x represents the scaled location of the impact ray which is 
deflected by a(x). The scaled deflection angle a is related to the true deflection angle a by 
a = D °^ ds a, where £ is an arbitrary scale factor. If the lens is symmetric around the line 
of sight, equation 7 reduces to a function of only the radial component x, and a lies along 
the radial direction x. 



2 f x 

a(x) = - xK[x')dx (8) 
x Jo 

Thus, to determine the lensing contribution of the cluster, we must first model n(x). 

From equation 3, we find that the three dimensional mass density at a radius r can be 
expressed as, 



p (r) mT (3 + (r/rc)2) (9) 

p9rav{ ) 47r/im p Gr c 2 (1 + (r/r c ) 2 ) 2 [ ) 

This quantity can be projected on the plane of the sky (normal to the line of sight), 
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resulting in the following expression for the surface mass density, £(#), as a function of the 
cylindrical radius x: 



3kT(3 {r 2 ^(r 2 + x 2 )(R 2 -x 2 ) + (2r 4 c + 2r 2 R 2 + r 2 x 2 + x 2 R 2 )<5>) 
^ = 2fxm p irG (r 2 + x 2 f/ 2 (r 2 c + R 2 ) (10) 

where, 

I R 2 - x 2 

$ = tan , 
and R is the radial extent of the cluster. 



\ r 2 + 



x 2 



Notice that the surface mass density of the cluster and the critical mass density both 
depend linearly upon H , thus implying that k does not depend on H . However, in our 
derivation of the convergence we infer the temperature of the cluster through an empirical 
luminosity - temperature relation and therefore our results do depend on the assumed value 
for /Jo- 
in determining k(x) we have examined a variety of models in which we used a 
combination of different cluster core parameters. We chose cluster limits ranging between 
120" and 280". We also varied the core radius between the values of 9 C = 5" and 9 C = 45" 
and use Hq values of 50 and 75 km s" 1 Mpc" 1 . Note that at the smaller radii, the model 
approaches the limit of a singular isothermal sphere. We find that the average value of k, 
at the location of the images, varies in the range of 0.07 to 0.21 as shown in Figure 6. The 
upper and lower bounds on k place respective upper and lower bounds on H in accordance 
with the models of Falco et al. (1991) and GN. For the estimation of H we have taken a 
self consistent iterative approach. We begin with an initial trial value for H and derive 
the total number of counts and cluster luminosity for an assumed core radius and cluster 
extent. Using the Lx - T relation, we derive the cluster temperature and from equations 
3 and 5 we obtain the total mass and convergence of the cluster, respectively. The GN 
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result (equation 4) provides an estimate for H . We reuse this value for H to initialize the 
iterative procedure and repeat it until the value for H converges (AH < 1). The result is 
very insensitive to the initial guess for H . The error for H is calculated by applying the 
self consistent iterative approach for cluster core radii ranging between 5" and 45", cluster 
radii ranging between 120" and 280" and changing the mean HRI count rate by ± 2(7. 
Assuming q = 0, we obtain 67 < H < 82 km s" 1 Mpc -1 for a time delay of 1.1 years. 

The effect of the cluster involves a rotation and an elongation of the separation vector 
connecting the two images. We may express the deflection angles for images A and B as 

— * — * 

O-a = k-aOoa and as = kb@ob, respectively, where the subscript on k indicates the point of 

— * — * 

evaluation and where do a and 6ob are the angular distances from the cluster center to the 

image locations of A and B, respectively. The calculated values for the deflections a a and 

a B are shown in Figure 7. The fractional elongation of the angular separation of images A 

I ab \ — \a^b' I 

and B produced by the cluster potential alone is - — — - and the calculated values for 

— * — * 

the percent elongation are presented in Figure 8; here AB and A'B' denote, respectively, 
the vector separations of the images with and without the cluster contribution. 

To calculate the convergences, ka and Kb, we use equation 8 with angular distances 
between images A and B and the cluster center given by xa = 26.57" and xb = 21.04", 
respectively, based on the published position (-13.7", -19.6") (Angonin-Willaime et al. 
1994) for the estimated cluster center of mass with respect to image B and the measured 
position of image A with respect to image B of (-1.25252" ,6.04662") (Gorenstein et al. 
1984). Again, we use a grid of models in which we vary the cluster core radius from 5" to 
45" and the cluster extent from 120" to 280". 



3. SPECTRAL ANALYSIS 



The quasar Q0957+561 also was observed with the ROSAT PSPC on 1991 November 
14 for 19006 sec, on 1993 November 8 for 4972 sec, and with the ASCA GIS on 1993 
May 9 for 29000 sec. A summary of information related to each of these observations is 
presented in Table 1. The spectral analysis for the 1991 November ROSAT observation was 
presented in Chartas et al. (1995). The 1993 November PSPC observation was centered on 
Q0954+556 leaving the quasar images located at an off-axis angle of about 47 arcmin. For 
the analysis we extracted the source counts within a circular region centered on Q0957+561 
with a radius of 6 arcmin. Background counts were extracted from several other circles 
of 6 arcmin radius with centers located at an off-axis angle of 47 arcmin, equal to that 
of Q0957+561. To model the spectral response of the PSPC, we used the September 
1994 PSPC detector response matrix pspcb_gain2_256.rmf provided by the ROSAT guest 
observer facility (GOF). The spectrum was grouped to obtain a minimum of 20 counts in 
each energy bin, allowing use of x 2 statistics. 

A simple redshifted power law plus galactic absorption due to cold material with solar 
abundances, provides an acceptable x 2 per degree of freedom (see Table 2). However, the 
resulting power law is slightly steeper (Aa u ~ 0.3) than expected for a radio loud quasar at 
a redshift of 1.4 (Wilkes & Elvis 1987; Shastri et al. 1993, Laor et al. 1997, Schartel et al. 
1996). 

In general the results obtained from the spectral fits to the 1993 November PSPC data 
are consistent with the results from the 1991 November PSPC observation of Q0956+561. 
In particular we find in both cases that an absorbed redshifted power law of photon index 
a v plus a thermal Raymond-Smith model provides an acceptable x 2 fit to the data and 
yields values for the fitted parameters a u , T e and N H that are close to what is expected for 
a radio loud quasar at z=1.4. The X-ray flux of the thermal component derived from the fit 
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to the PSPC spectrum of Q0957+561 is about 5 x 1CT 13 ergs s _1 cm -2 , at least a factor of 
five greater than the X-ray flux of the cluster as derived from the ROSAT HRI observations 
(see Table 3). This result implies that the observed ROSAT PSPC thermal component 
most likely originates from a source other than the cluster (see, e.g., Chartas et al. 1995). 

To constrain the model parameters even further we performed simultaneous fits to 
the 1991 November and 1993 November PSPC data. We allowed the normalization of the 
power law and thermal component of each observation to vary independently. 

A comparison of the 0.2 - 2keV X-ray fluxes between the 1991 and 1993 ROSAT PSPC 
observations of Q0957+561 implies a decrease in flux over this interval by (25 ±9)% with 
the quoted error being at the 90% confidence level. 

Q0957+561 was also observed with the ASCA satellite during a pointed observation of 
NGC 3079 on 1993 May 9 for 29072 sec. In this observation Q0957+561 is located 14.1' off 
axis in the GIS2 image and 9.28' off axis in the GIS3 image just within the field of view of 
the GIS2 and GIS3, but unfortunately outside the solid-state imaging spectrometer's (SIS's) 
field of view. The source counts were extracted from circles of radius 3.5 arcmin centered 
on Q0957+561. The background was extracted from several other circles of the same radius 
centered at an off- axis angle equal to that of Q0957+561, but far away in azimuth. For 
the data analysis of the GIS2 and GIS3 spectra we used the spectral response functions 
gis2v4_0.rmf and gis3v4_0.rmf, respectively, provided by the ASCA GOF. The spectral data 
were grouped to obtain a minimum of 20 counts per bin. The observed ROSAT PSPC 
and ASCA GIS spectra are shown in Figure 9. The few source counts obtained in the 
ASCA observation prevent us from fitting any complex models; however, a comparison of 
the spectral slope of a simple power law model fitted to the 0.5-6 keV GIS data to that of 
the slope obtained from fits to the 0.2-2 keV ROSAT data was possible. In particular a 
simultaneous fit of a simple absorbed power law model between 0.5 and 5keV to the GIS2 
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and GIS3 spectra of Q0957+561 yields a photon index a u = 1.6 ± 0.1, quite typical for 
radio loud quasars at z—lA; a simultaneous fit between 0.2 and 2keV of a similar model 
to the November 1991 and November 1993 ROSAT PSPC spectra gave a photon index of 
2.27 ± 0.03. This difference may indicate a break in the spectrum at some energy above 2 
keV. Alternatively, including a thermal component (fits 13 and 14) results in values of the 
photon index of the power law component consistent with that for radio loud quasars. 

4. DISCUSSION and CONCLUSIONS 

The need for a detailed measurement of the mass distribution of the cluster lens has 
been emphasized by Falco et al.(1991) and by GN (1996) in analyses of lens models for 
Q0957+561 that incorporate the cluster contribution through a convergence parameter K,(r). 
The most important result of the present analysis of the ROSAT HRI data of Q0957+561, 
that we hope will eventually lead to an accurate determination of the mass distribution of 
the cluster lens, is the 3a detection of X-rays from the cluster of galaxies that contains the 
principal lens galaxy Gl. 

As indicated in Figure 3, the uncertainty in our estimated values for the mass 
distribution of the cluster spans a large range as do the consequent estimates of the 
convergence parameter of the lens (Figure 6). The main contributors to this uncertainty 
are the lack of knowledge of the shape of the cluster, the core radius of the cluster, the 
temperature profile of the hot gas, the cluster extent, and the location of the cluster 
center as well as the sensitivity to the value of H used for deriving the cluster hot gas 
temperature from the X-ray luminosity. The present analysis yields a total cluster mass 
in the range 1.5 - 3.2 x 1O 14 M0, and does not include allowances for errors in several of 
the above characteristics, but does consider the effects of cluster core radii between 5" and 
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45", and cluster extents between 120" and 280". The calculations of the X-ray luminosity 
were performed for Hubble constants of 50 and 75 km s^ 1 Mpc -1 . The average convergence 
parameter k = (ka + «b)/2 of the cluster at the image locations A and B was found to lie 
between 0.07 and 0.21. 

One of the other assumptions in the present calculation is that the temperature of the 
cluster follows the luminosity temperature relation of David et al. (1993) for a large sample 
of clusters of galaxies. The derived cluster temperature depends weakly on the cluster 
X-ray luminosity. An error of a factor of 2 in the X-ray luminosity derived from the HRI 
data will propagate as an error of about 30% in the cluster temperature and about 7% (for 
k = 0.2) in the derived value of the Hubble constant. The residual rms scatter in the L x - 
T relation along the temperature axis is <Ji og T = 0.104/2.1 according to Markevitch (1998) 
and propagates as an error of about 3% in the derived value of the Hubble constant. 

Our calculations for the deflection angles produced from considering the influence of 
the cluster alone (not including Gl) found the deflection for images A and B to lie in the 
range of 2" to 16" as shown in Figure 6; the differences in these deflections for images A 
and B are much smaller. The resulting percent contribution, ASab, to the observed 6."1 
separation of images A and B due to the cluster alone, lies in the range of 4% to 22%. This 
range includes only the effects of varying the cluster core radius between 5" and 45" and the 
cluster limit between 120" and 280". This result is in contrast with the commonly accepted 
notion that the large separation of the two images of Q0957+561 is an indicator that the 
lens contains a cluster of galaxies. To test the sensitivity of this result to the adopted cluster 
center position, we computed the cluster contribution to the image separation for cluster 
center locations within a distance of 5" of the location derived from the optical observations 
by Angonin-Willaime et al. (1994). As expected, when the cluster center is moved closer to 
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the locations of images A and B, the derived cluster contribution to the image separation 
becomes more significant reaching a value as high as ASab — 31% for a cluster center at a 
distance of about 16" from image B. 

Recently the mass distribution of the lens of Q0957+561 was determined by Fischer 
et al. (1997) based on the distortion of background galaxies produced by the lens. Their 
estimate of 3.7 ± 1.2 x 1O 14 M0 within a radius of IMpc is consistent with the X-ray 
estimates. As pointed out by Schneider & Seitz (1995) however, the Kaiser & Squires (1993) 
mass reconstruction technique also is insensitive to constant density sheets of matter. 

An interesting conclusion from the mass reconstruction analysis by Fischer et al. 
(1997) is the relatively low value of 5" derived for the cluster core radius. The cluster core 
radius, r cF , as defined in Fischer et al. (1997), however, differs from the cluster core radius, 
r c , used in this paper and defined in equation 2. The relation between r c and i c p is given 
approximately by r cF = -0.987 + 0.685r c for 5" < r c < 50". We will have to await future 
X-ray observations of Q0957+561 to compare more accurately X-ray and weak- lens results 
for the cluster core radius. 

The present ROSAT detection of cluster emission suggests that future X-ray 
observations will greatly add to our understanding of Q0957+561. AXAF observations 
should provide direct determination of the temperature profile, the core radius and the 
shape of the cluster lens, reveal possible dumpiness in the total mass profile, and ultimately 
provide a tighter constraint on the Hubble constant. 

We are grateful to Emilio Falco for providing many useful comments and to Donald 
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Fig. I. — Composite iso-intensity plot of three co-aligned ROSAT HRI observations of 
Q0957+561. Contour levels are expressed as a percentage of the peak value. 

Fig. 2. — The total counts from the ROSAT HRI, N C i, within a circular region centered on 
Q0957+561 and originating from cluster emission. This number was estimated by extracting 
from the HRI image of the Q0957+561 system the number of counts within a selected annulus 
centered on Q0957+561 and applying a correction for the expected counts residing outside 
this annulus. For this plot, we used a cluster core radius of 45" and a cluster extent of 280". 

Fig. 3. — Estimated total mass of the cluster containing the Gl galaxy within a radius of 
IMpc over a range of typical cluster core radii, for q = and for H values of 50 and 75 
km s" 1 Mpc -1 . 

Fig. 4. — Conceptual diagram of the cluster's contribution to the deflection effect. Images A 
and B are the observed lensed images of Q0957+561 while images A' and B' are the expected 
images for deflection produced from the Gl galaxy alone. 

Fig. 5. — Basic lens geometry of the system Q0957+561. The shown deflection angles c?a 
and as are due to the cluster contribution alone. 

Fig. 6. — The local convergence parameter k, averaged over locations of images A and B, 
as estimated for cluster core radii ranging between 5" and 45", and cluster limits ranging 
between 120" and 280". The calculations of the X-ray luminosity were performed for Hubble 
constants of 50 and 75 km s _1 Mpc -1 and for q = 0. This range of uncertainty in 
the convergence parameter k does not include possible systematic errors arising from the 
derivation of the cluster temperature through the use of the luminosity-temperature relation 
for clusters of galaxies or from other assumptions, such as the shape of the cluster. 

Fig. 7. — Calculated values for the deflection angles a a and «b produced by the cluster. 
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The cluster core radius, the cluster extent, and H were varied within the ranges specified 
in the caption for Figure 6. For other relevant assumptions see text. 

Fig. 8. — The percent elongation of the angular separation of images A and B produced by 
the cluster potential alone. The cluster core radius, the cluster extent and H were varied 
within the ranges specified in the caption for Figure 6. For other relevant assumptions see 
text. 

Fig. 9. — Simultaneous spectral fit to ROSAT PSPC data taken on 1991 November 15 and 
1993 November 8, and to ASCA GIS data taken on 1993 May 9. The model used for this 
fit incorporates a redshifted power law plus a thermal Raymond-Smith model and galactic 
absorption (see text). 
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TABLE 1 

ROSAT HRI, PSPC and ASCA GIS Observations of Q0957+561 



Observation 


Instrument 


Observation 


Exposure 


Off axis 


Net Counts 


Date 




File name Id 


sec 


angle 
arcsec 




1991 May 1 


ROSAT HRI 


rh700100 


20847 


45 


473 ± 30 


1992 October 20 


ROSAT HRI 


wh600411n00 


20676 


740 


812 ± 34 


1992 October 22 


ROSAT HRI 


rh700889n00 


18548 


50 


800 ± 33 


1994 April 11 


ROSAT HRI 


rh701294n00 


25083 


50 


1077 ± 39 


1991 November 15 


ROSAT PSPC 


rp300319 


19006 


993 


3677 ± 70 


1993 November 8 


ROSAT PSPC 


rp701401n00 


4472 


2821 


449 ± 27 


1993 May 9 


ASCA GIS2 


ad60000000g200470 


29072 


846 


140 ± 20 


1993 May 9 


ASCA GIS3 


ad60000000g300470 


29072 


557 


121 ± 18 



NOTE- A description of the format for the ROSAT high resolution instrument (HRI), ROSAT position sensitive proportional 
counter (PSPC) and ASCA gas imaging spectrometer (GIS) file names is available at the ROSAT and ASCA guest observer 
facility web pages. 
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TABLE 2 

Model Parameters Determined from Spectral Fits to the 
ROSAT PSPC and ASCA GIS Spectra of Q0957-561 



Fit 


Model" 


Instruments 6 






7V H (z = 0) 
10 20 cm- 2 


kT e (keV) 


X 2 /dof(dof) 


1 


1 


PSPC t2 




2 36+ 011 


0.87 




0.95(13) 


2 


2 


PSPCt2 




2 g7 +0.45 


2 16+°" 




0.80(12) 


3 


3 


PSPCt2 




, 71 +0.44 
L - 1 1 -0.7S 


0.87 


43+0.11 


0.73(11) 


4 


1 


PSPC t 2 + PSPCti 




9 9 7 +0.028 
z - z -0.027 


0.87 




0.88(63) 


5 


2 


PSPC t 2 + PSPCti 




z - oo -0.08 


1 44+0-22 




0.69(62) 


6 


3 


PSPC t 2 + PSPCti 




2 06+ - 07 


0.87 


n 45+0-08 


0.62(60) 


7 


4 


PSPC t2 + PSPCti 




2 14+ - 85 


99+ - 44 


42+ ' 09 
-0.28 


0.63(59) 


8 


1 


GIS2 




i oq+0-12 
-"-•^ -0.022 


0.87 




1.85(12) 


9 


1 


GIS3 




1 m+0.19 


0.87 




1.34(7) 


10 


1 


GIS2 + GIS3 




i.oo_ 12 


0.87 




1.99(20) 


11 


1 


PSPC t 2 + PSPCti 
GIS2 + GIS3 


+ 


9 9CT+0.07 
z - ZJ -0.08 


0.87 




1.42(84) 


12 


2 


PSPC t 2 + PSPCti 
GIS2 + GIS3 


+ 


9 oo+0.07 


1 04+ - 15 




1.42(83) 


13 


3 


PSPC t 2 + PSPCti 


+ 














GIS2 + GIS3 




1 89+ - 08 


0.87 


n 41+0-05 


1.09(79) 


14 


4 


PSPC t 2 + PSPCti + 














GIS2 + GIS3 




1 58+ ' 14 


42+ - 14 

U - 4Z -0.14 


u -°-o.i 


1.01(78) 



NOTES- 

a Model 1 incorporates a redshifted power law plus absorption due to cold material at solar abundances fixed to the galactic 
value. Model 2 incorporates a redshifted power law plus absorption due to cold material at solar abundances set as a free 
parameter in the fit. Model 3 incorporates a redshifted power law plus thermal Raymond-Smith model plus galactic absorption. 
Model 4 incorporates a redshifted power law plus thermal Raymond-Smith model plus absorption due to cold material at solar 
abundances set as a free parameter in the fit. 

6 PSPCti and PSPC t 2 correspond to the ROSAT PSPC observations made on 1991 November 15 and on 1993 November 8, 
respectively. 
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TABLE 3 

X-ray Fluxes and Luminosities of the Thermal and Power Law Components 



Fit 


Instrument 


fa 

J X , Thermal 

10~ 12 erg cm~ 2 s _1 


fa 

J X ,P ower 

10~ 12 erg cm -2 s _1 


^X,Thermal 

10 46 erg s- 1 


^ X , Power 

10 46 erg s- 1 


6 


PSPCn 


0.52 ±0.3 


1.03 ±0.19 


1.29 


1.68 


6 


PSPCt2 


0.57 ±0.3 


0.68 ±0.19 


1.42 


1.11 



NOTES- 

a Absorbed flux between 0.2 - 2 keV. 

b Luminosity between 0.2 - 2 keV. No corrections due to magnification from lensing have been made. 
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